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ABSTRACT 

Timing observations of 40 mostly young pulsars using the ATNF Parkes radio telescope 
between 1990 January and 1998 December are reported. In total, 20 previously unre- 
ported glitches and ten other glitches were detected in 11 pulsars. These included 12 
glitches in PSR J1341— 6220, corresponding to a glitch rate of 1.5 glitches per year. We 
also detected the largest known glitch, in PSR J1614-5047, with Av g /v « 6.5 x 10~ 6 
where v = 1/P is the pulse frequency. Glitch parameters were determined both by 
extrapolating timing solutions to inter-glitch intervals and by phase-coherent timing 
fits across the glitch(es). These fits also gave improved positions and dispersion mea- 
sures for many of the pulsars. Analysis of glitch parameters, both from this work and 
from previously published results, shows that most glitches have a fractional amplitude 
Aug/v of between 10~ 8 and 10 -6 . There is no consistent relationship between glitch 
amplitude and the time since the previous glitch or the time to the following glitch, 
cither for the ensemble or for individual pulsars. As previously recognised, the largest 
glitch activity is seen in pulsars with ages of order 10 4 years, but for about 30 per 
cent of such pulsars, no glitches were detected in the 8-year data span. There is some 
evidence for a new type of timing irregularity in which there is a significant increase in 
pulse frequency over a few days, accompanied by a decrease in the magnitude of the 
slowdown rate. Fits of an exponential recovery to post-glitch data show that for most 
older pulsars, only a small fraction of the glitch decays. In some younger pulsars, a 
large fraction of the glitch decays, but in others, there is very little decay. Apart from 
the Crab pulsar, there is no clear dependence of recovery timescale on pulsar age. 
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1 INTRODUCTION 

Timing observations have shown that pulsars are remarkably 
stable clocks (e.g. Kaspi, Taylor & Ryba 1994). However, 
they are not always predictable. Two types of instability 
have been observed. In the first, the pulsar period, P, or 
frequency v = 1/P, varies in an apparently random fash- 
ion, fluctuating on timescales of days, weeks, months and 
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years ( Cordes fc Downs 1985 ). The intensity of the fluctua- 
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tions can be described by an 'activity parameter' based on 
the amplitude of the second frequency-derivative term in a 
Taylor-series fit to a set of pulse times of arrival (TOAs): 



(1) 



A(t) = log ( ^-r 



where t is the length of the data span in seconds ( Axzouma- 



nian et al. 1994). Conventionally t = 10 s is adopted, giving 



the parameter As . In most young pulsars, the magnitude of 
v is far in excess of that expected from secular slowdown. 
Furthermore , timing noise is normally very 'red' (Cordes & 



Downs 1985 ) and so the second derivative term is a good in- 
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dicator of the level of noise. Observations of large samples of 198£ ; Alpar et al 



pulsars flCordcw fc Downs 1985 ; D'Amico ct al. 199£ ) show 
that the activity parameter is positively correlated with the 
value of P, the first time-derivative of the pulsar period. 

In the second type of period instability, known as a 
'glitch', the pulse frequency has a sudden increase which 
typically has a fractional amplitude Av g jv in the range 10 -8 
- 10~ 6 . These glitches are unpredictable but typically occur 
at intervals of a few years in young pulsars. Coincident with 



1995 



Ruderman, Zhu & Chen 1998) are 



the glitch, there is often an increase in the magnitude of 



the frequency derivative, typically by ~ 1 per cent, which 
sometimes decays with a timescale of weeks to years. The 
time- dependence of the pulsar frequency after a glitch is 
generally well described by the following function: 



generally consistent with a rather small superfluid fraction, 
I s /I ~ 1CP 2 , participating in the glitch activity. In the Ru- 
derman et al. (1998) model, permanent changes in v may 
result from a change in the magnetic field configuration as- 
sociated with the crustal cracking at the time of a glitch. 

Until the last few years, the number of known glitches 
and glitching pulsars was modest. However, with re latively 
high-frequency searches at low Galactic latitudes ( Clifton 
Johnston et al. 1992] ; Kaspi et al. 1992) discov 



et al. 1992 



v{t) = u {t) + Av g [l - Q(l - exp(-t/r d )] + Av v t 



(2) 



where vo(t) is the value of v extrapolated from before the 
glitch, Av g = Ai>d + Av v is the total frequency change at 
the time of the glitch (t = 0), where Avd is the part of the 
change which decays exponentially and Av v is the perma- 
nent change in pulse frequency, Q — Aud/Au g , Td is the 
decay time constant and Ai> v is the permanent change in z> 
at the time of the glitch. The increment in frequency deriva- 
tive is given by 



Av(t) = Av d exp(-t/r d )~ 



Ai/p = 



~ QA " 3 exp(-t/T d )+Au p . 

Td 



The middle term of equation |^ was originally derived 
on the basis of the two-component model for glitch recovery 
(Baym et al. 1969[ ), where the initial frequency jump was 



due to a 'starquake' or sudden change in moment of inertia 
of the solid crust. In this model, the parameter Q — I s /I c 
where I s and I c are the moments of inertia of the superfluid 
component and the crust, respectively. 

In later versions of the theory, the frequency jump is 
due to a su dden unpinning of vortex lines in crustal neutron 
superfluid ( Anderson fc Itoh 1975 ). Post-glitch relaxation 
may be due to 'vortex creep', that is, a slow drift of the 
vortices across the crustal lattice (e.g. Alpar, Cheng & Pines 
1989) or drift of t he crustal lattic e itself with the vortices 
rem aini ng pinned (Eludcrman 1991; Ruderman, Zhu & Chen 
199s£] 

Depending on the internal temperature of the star and 
other f actors, the response to a glitc h may be linear or non- 
linear (Alpar, Cheng & Pines 1989). In the linear regime, 
the frequency jump at t = decays exponentially with time 
constant Td. The fractional change in frequency derivative 
at the jump is 



Av a jv = (I s /I)(Au s /T d ), 



(4) 



where Auj s is the change at the time of the glitch in the 
rotational lag of the superfluid, and uj s — 2n{v s — v), where 
v s is the rotation frequency of the neutron superfluid. For 
glitches which are large compared to the steady-state lag, 
Alu 3 ~ — 2tvAv and the change in spin-down rate can be 
large compared to I a /I- 

In the non-linear regime, there is essentially no relax- 
ation (Q » 0) and 



AVgjv = I a /I, 



(5) 



giving a permanent change Az> p in v (Eqn |3j). Analyses of 
data for the Vela pulsar and other pulsars (Alpar et al 



ering a much larger sample of young pulsars, the number has 
increased significantly, with 21 pulsars having 46 glitches be- 
ing listed in the recent review by Lyne (1996). In this paper 
we describe timing observations of a sample of 40 mostly 
young pulsars using the Parkes radio telescope which show 
a total of 30 glitches in 11 pulsars. 



2 OBSERVATIONS AND ANALYSIS 

Observations were made using the Parkes 64-m radio tele- 
scope between 1990 January and 1998 December at frequen- 
cies around 430, 660, 1400 and 1650 MHz. The frequencies 
below 1 GHz were used only occasionally to improve the de- 
termination of dispersion measures (DMs). At all frequen- 
(3) cies, cryogenic dual-channel receivers were used. Two sepa- 
rate back-end systems were used to provide the frequency 
resolution necessary for dedispersing. The early observations 
used filterbanks and a one-bit digitizer system constructed 
at Jodrell Bank; a 2 x 256 x 0.125 MHz system was used 
at 430 MHz, a 2 x 128 x 0.25 MHz system was used at 660 
MHz and at the higher frequencies a 2 x 64 x 5 MHz sys- 
tem was used. Further details of the filterbank data acquisi- 
tion system may be found in Manchester et al. (1996). From 
1994 July, data were acquired using a correlator system con- 
structed at Caltech (Navarro 1994). This system used two- 
bit digitization and an autocorrelator with 2 x 256 lags over a 
maximum bandwidth of 128 MHz. From mid-1995, the lags 
were split between two frequency bands which gave simul- 
taneo us observations at r adio frequencies of 1400 and 1650 
MHz ( jSandhu ct al. 1997j ). 

For the filterbank systems, data were folded syn- 
chronously with the pulsar period using off-line programs 
to give mean total intensity pulse profiles. In the correlator 
system this folding was performed in a hardware integrator- 
having 1024 bins across the pulsar period. For both systems, 
integration start times were established to better than \ps 
using time signals from the Observatory clock system. This 
was related to UTC(NIST) using a radio link to the NASA 
DSN station at Tidbinbilla and clock offsets kindly provided 
by the Jet Propulsion Laboratory, Pasadena. In subsequent 
analysis, the data were summed to form 8 or 16 frequency 
subbands and time subintegrations of 60 or 90 s duration and 
stored on disk. These files were then summed in frequency 
and time using the best available period and DM informa- 
tion to form a single profile for each observation. This was 
then cross-correlated with a standard profile for the pulsar 
to give a pulse time of arrival (TOA). 

A total of 40 pulsars were monitored during the pro- 
gram. The J2000 and B1950 names of these pulsars, their pe- 
riods and characteristic ages, and the dates spanned by the 
observations are listed in Table 0. Depending on the strength 
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of the pulsar, observation times for each TOA were between 
2 and 12 min. Observations were obtained at intervals of be- 
tween 2 and 6 weeks for most of the pulsars in Table ^, with 
the known glitching pulsars being observed more frequently. 
The TO As resulting from this program were analysed using 
version 11.3 of TEMPO0 which includes provision for fit- 
ting the glitch parameter s given in Equa tion ^. The DE200 
solar-system ephemeris (|Standish 1990) was used to con- 



vert TOAs to the solar-system barycentre. The sixth and 
seventh columns of Table ^ give the rms value of the TOA 
uncertainty (dependent on the pulse period, pulse shape and 
signal/noise ratio) and the rms residual after fitting for just 
pulse frequency and its first derivative. For glitching pulsars, 
the largest glitch- free interval was used for this fit. These 
residuals are dominated by the effects of timing noise, and 
so they are an indication of its amplitude. An exception to 
this is PSR J1513-5908, which has significant frequen cy sec- 
ond an d th ird derivatives due to secular slowdown (Kaspi 
et al. 1994); in this case, the fit giving <tt included the first 
three frequency derivatives. 

For several pulsars, the position and DM were mea- 
sured with comparable or better precision than was previ- 
ously available. Positions were obtained using the so-called 
'pre- whitening' method (Kaspi et al. 1994). In this method, 
the position is fitted along with sufficiently many frequency- 
derivative terms and in some cases, glitch terms, to 'absorb' 
the timing noise and give an approximately Gaussian distri- 
bution of timing residuals. DMs were measured from short 
sections of data where there were multi-frequency observa- 
tions. This ensures that contamination by long-term timing 
noise is not a problem and that error estimates are realistic. 

For glitches where the amplitude of the glitch and the 
interval between the last pre-glitch observation and the first 
post-glitch observation are not too large, the epoch of the 
glitch can be determined by requiring that the pulse phase 
be continuous over the glitch. In this situation, TEMPO 
gives an estimate of the glitch parameters, including the 
epoch, and their uncertainty. Where post-glitch decay pa- 
rameters are estimated, these refer to the long-term decay, 
typically with timescales of hundreds of days. The resolu- 
tion of our observations is generally not sufficient to detect 
the more rapid post-glitch recoveri es observed in s o me pu 



sars, fo r examp le, the Vela pulsar (Flanagan 1990 



al. 199C) 



McCul- 



loch et 

For larger glitches where there may be one or more turns 
of phase in the residuals between the bounding observations, 
the glitch epoch is uncertain. Other glitch parameters are af- 
fected by this uncertainty. To estimate uncertainties in this 
case, the following procedure was adopted. The glitch epoch 
was taken to be halfway between the bounding observations, 
with an uncertainty of half their separation. The increments 
Av g and Av g were then computed for an assumed glitch 
epoch t g by extrapolating the pre- and post-glitch fits to t g 
and taking differences. Uncertainties were similarly extrap- 
olated and quadrature sums taken for the uncertainties in 
the increments. This was done separately for t g at the two 
bounding epochs. The final error estimates were then the 
quadrature sum of the difference between the increments at 



either end of the data gap and the larger of the uncertainties 
in the increments. 



3 RESULTS 

Table ^| lists parameters for the 11 pulsars for which glitches 
were detected. The pulsar J2000 positions and DMs used in 
or derived from the analysis are given. Uncertainties in the 
last digit quoted are given in parentheses. These and other 
quoted uncertainties from TEMPO fits are twice the formal 
rms error. Position fits were generally to the largest avail- 
able data span not too strongly affected by post-glitch post- 
glitch recovery. As described above, the data spans used to 
determine the positions and their errors were 'pre-whitened' 
to eliminate the effects of period noise from the positions 
and their errors. The fifth and sixth columns give the data 
span used when fitting the position and the final rms timing 
residual. References for the position and DM are given in 
the final column. 

A total of 30 glitches were observed in these 11 pul- 
sars. Independent fits to the pre-glitch, inter-glitch and post- 
glitch timing data are given in Table |j| Except for short sec- 
tions of data, the fits include a v term. In every case, this 
term is dominated by recovery from a previous glitch; the v 
from the long-term secular slow down is negligible by com- 
parison. Higher-order frequency derivative terms were not 
fitted, so the rms residuals given in right-most column re- 
flect the presence of random period irregularities, especially 
for longer data spans. 

Glitch parameters are listed in Table ^[ If marked by *, 
the glitch epoch is determined by requiring phase continuity 
across the glitch. The next two columns give the fractional 
steps in v and v at the glitch determined by extrapolating 
the pre- and post-glitch solutions (Table w to the glitch 
epoch, with uncertainties determined as described in Section 
2. Parameters in the remaining columns were determined 
using TEMPO to fit phases across the glitch. 

Ten of these glitches have bee n previously rep orted: 
single glitche s in PSR J0835-4510 ([Flanagan 1996|), PSR 
J1709-4428 (|johnston et al. 1995^ |Shemar fc Lyne 1996D , 
PSR J1731-4744 (p'Alessandro fc McCulloch 1997^ ~and 
PSR J1801-2451 ( I Lyne et al. 1996a| ) and three glitches in 
PSR J1341-6220 ( |rCaspi et al. 1992|; |5hemar fc Lyne 199fj) 



and PSR J1801-2306 ( [Kaspi et al. 1993[ [Shcmar fc Lyne 



199f). We reanalyse these glitches for completeness and con- 



sistency with the results for the other pulsars. 

In the following sections we discuss each pulsar in turn. 



3.1 PSR J0835 4510, PSR B0833 45, the Vela 
pulsar 

The Vela pulsar is well known to suffer many giant glitches 



(|Cordcs, Downs fc Krausc-Polstorff 1988 



VIcCulloch et al 



t See http://pulsar.princcton.edu/tcmpo/ for a description of 
TEMPO. 



199C; Lyne et al. 1996b) and is being regularly monitored 
at a number of observatories. At Parkes, we are unable to 
make such regular and frequent observations. In this paper, 
we describe observatio ns of the latest glitch which has only 
been briefly reported ( Flanagan 1996 ). 

Table H gives the results of fitting for v and its first 
two time derivatives, both before and after the glitch. The 
pre-glitch fit was obtained from an approximately 1-year 
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Table 1. Pulsars monitored for glitch activity 
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data sp an and shows a significant v resulting from previous haviour is sim ilar to that seen in previous Vela glitches (Lync 



glitches. The post-glitch fit given in Table 3 is for an 82-day 
span commencing 18 days after the glitch. A timing model 
with polynomial terms up to v is not a good fit to longer 
data spans. 

Fig. |l] shows the time dependence of the frequency resid- 
ual Av and of v around the time of the glitch. The Au values 
plotted are differences between the values of v obtained from 
independent fits to short sections of data, typically of span 
20 - 30 days, and those determined from the predictions 
of the pre-glitch model given in Table |^. This plot shows 
an approximately exponential recovery in v after the glitch, 
followed by an approximately linear increase in v. This be- 



et al. 1996b) 



Flanagan (1997) gives Auju = 2.15(2) x 10~ 6 and an 
epoch for the glitch of 1996, October 13.394 UT, correspond- 
ing to MJD 50369.394. Fractional changes in u and i> at the 
time of the glitch, obtained by extrapolating the pre-glitch 
and post-glitch fits (Table ^) to the glitch epoch are given 
in the fifth and sixth columns of Table ^. Results of fitting 
the exponential model (Equation |^) to the interval from 350 
days before the glitch to 200 days after are also given in 
Table [l|. In this fit, p, u and v were held at their pre-glitch 
values (Table |^) . The small rms residual shows that the ex- 
ponential model with about 40 per cent of the initial glitch 
in frequency decaying on a timescale of about 900 d is a 
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Table 2. Parameters for glitching pulsars 
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18:01:19.829(9) 


-23 


04:44 


2(2) 


1801- 


-2451 


18:01:00.223(7) 


-24 


51:27 


1(1) 


1803- 


-2137 


18:03:51.35(3) 


-21 


37:07 


2(5) 



DM 


MJD Range 


aw 


Ref. 


(cm -3 pc) 






(lis) 


Posn,DM 


68.094(4) 






— 


1,5 


129.10(1) 


49043- 


-50536 


1150 


5,1 


271.01(4) 


49176- 


-50402 


1080 


5.2 


223.26(3) 


49708- 


-51155 


785 


5.3 


717.3(6) 


48874- 


-49888 




5,5 


582.8(3) 


50269- 


-50778 


736 


5,1 


75.69(5) 


48928- 


-51156 


211 


5,1 


123.33(2) 


49415- 


-50704 


351 


5,5 


1074(6) 








1,1 


289(1) 








1,1 


233.9(3) 








1,1 



a Reference: 1. Taylor, Manchester & Lyne (1993) ; 2. Kaspi et al. (1996) ; 3. D'Amico et al. (1997) ; 4. Frail, 
Kulkarni & Vasisht (1993) ; 5. This paper. 

b Position is for epoch MJD 41380 with proper motion [i a = —48 mas yr — 1 , [ig = 35 mas yr~ 1 
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(a) 



o 



360 /-is and a systematic quartic term in the residuals. For 
the exponential fit, most of the rms residual comes from the 
first few post-glitch points which are positive, indicating an 
additional short-term recovery. However, it was not possi- 
ble to fit for the parameters of this. The fitted glitch epoch, 
given in the third column of Table ^, is slightly earlier than 
that given by Flanagan (1997); this is consistent with the 
presence of more rapid post-glitch relaxation than modelled. 



(b) 
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Figure 1. The 1996 October glitch of PSR J0835-4510: varia- 
tions of (a) frequency residual l\v relative to the pre-glitch so- 
lution, (b) an expanded plot of Au where the mean post-glitch 
value has been subtracted from the post-glitch data, and (c) the 
variations of v. The epoch of the glitch is indicated by a small 
arrow near the base of each plot. 



very good representation of the post-glitch relaxation, at 
least over the data span fitted. Fitting of v, v and v to the 
post-glitch data instead of the three parameters of the ex- 
ponential fit gives a substantially worse fit with rms residual 



3.2 PSR J1048 5832, PSR B1046-58 

This young pulsar was discovere d in a high-frequency s urvey 
of the southern Galactic plane (Johnston et al. 1992). Ev- 
idence for gamma-ray pulsations has been found by Kaspi 
et al. (2000). The pulsar suffered two large glitches within 
the observed data span, the first of size Avg/v ~ 3 x 10~ 6 
in 1993 February (MJD ~ 49034) and the second of size 
~ 0.7 x 10~ 6 in 1997 December (MJD ~ 50788). Observed 
frequency residuals are shown in Fig. ^. The expanded plot 
of Av (Panel b) shows a third and much smaller glitch about 
100 d before the large first glitch. Fits to the inter-glitch in- 
tervals (Table |§|) give evidence for large fluctuations in the 
period. Fig. |^ gives timing residuals for the data span fol- 
lowing the largest glitch, showing quasi-random fluctuations 
with a timescale of a few hundred days. 

Table ^ gives an improved position for PSR J1048— 5832 
derived from the data following the largest glitch. Most of 
the systematic oscillation was removed by fitting up to the 
twelfth pulse frequency derivative at the same time as the 
position was determined. A fit to the whole data set, in- 
cluding parameters for the glitches, gave the same position 
within the combined errors. This position differs by 8" from 
that given by Johnston et al. (1995); this difference can be 
attributed to the different sets and the presence of the strong 
period irregularities. Recently, Stappers et al. (1999) have 
used the Australia Telescope Compact Array to determine 
an interferometric position for the pulsar: R.A. (J2000) 10 h 
48 m 12604 ± 0:008, Dec. (J2000) -58° 32' 03'.'75 ± 0'.'05. 
This position agrees with the timing position in R.A., but 
differs by 2'.'45 ± 0'.'80 in declination, probably as a result 



5 



Table 3. Pulse frequency parameters for glitching pulsars 



PSR J Int. v 









(s- 1 ) 


0835- 


-4510 


-1 


11.1964839822(3) 






1- 


11.1955738808(10) 


1048- 


-5832 


-1 


8.087302320(3) 






1-2 


8.0869925648(4) 






2-3 


8.0865824929(8) 






3— 


8.08599240181(15) 


1105- 


-6107 


-1 


15.8248916396(18) 






1-2 


15.8246473365(3) 






2- 


15.8245521708(20) 


1123- 


-6259 


-1 


3.68413905495(20) 






1— 


3.68413613509(3) 


1341- 


-6220 


-1 


5.173922035(4) 






1-2 


5.1737636204(9) 






2-3 


5.1735752982(4) 






3-4 


5.1733890567(14) 






4-5 


5.1732030931(11) 






5-6 


5.173030116(14) 






6-7 


5.1729425780(5) 






7-8 


5.1728338894(5) 






8-9 


5.172758037(3) 






9-10 


5.1726390106(3) 






10-11 


5.1724946720(3) 






11-12 


5.1723883167(8) 






12- 


5.1722422270(5) 


1614- 


-5047 


-1 


4.317954008(15) 






-1 


4.3177762511(8) 






-1 


4.3175461311(4) 






-1 


4.317231719(4) 






_1 


4 3169526367('20 > ) 






1_ 


4 31 67356028r6"l 






1- 


4.3163251092(3) 






1- 


4.3160405830(14) 


1709- 


-4428 


-1 


9.7612713776(3) 






1- 


9.7608645977(8) 






1- 


9.7599781350(7) 


1731- 


-4744 


-1 


1.20510332177(3) 






1-2 


1.20508592522(8) 






2— 


1 205067861 94('6~) 


1801- 


-2304 


-1 


2.4051602072(9) 






1-2 


2.40512168588(17) 






2-3 


2.40507531889(14) 






3-4 


2.40505493478(18) 






4- 


2.4050229993(3) 


1801- 


-2451 


-1 


8.007176212(3) 






1-2 


8.0065388371(19) 






2- 


8.0060495257(7) 


1803- 


-2137 


-1 


7.4832982197(3) 






1-2 


7.4832434845(19) 






2- 


7.4831211339(4) 







Epoch 


(10" 12 s- 2 ) 


(10- 24 s" 3 ) 


(MJD) 


1 ^ ^81 Af\( A\ 
— 10.0014^4 ) 


O 1 O^IU ) 


OUIOO.UU 


— 10.01 ( ou^D ) 


Gfic: ( 1 o\ 
yoo^lz ) 


0Uo4 ( .UU 


-6.27373(7) 


72(10) 


48419.00 


ft 971 Q/A\ 

— o. z ( iy i ft ) 




/iQQQi nn 
4oyy i .uu 


-6.28460(3) 


146.7(14) 


49790.00 


-6.29765(7) 


390(40) 


50889.00 


Q Q^QI Q(^\ 
— O.uOO Luyi) J 


-4y^o; 


{ oy .uu 


Q QfiQI Q/l r;\ 




c;n^i r nn 

OUOID.UU 


— o.yoozo^iu ) 


— 04 lo ) 


c;n7QA nn 
ou / y4.uu 


C\ fi71 OC/cN 

— U.U J 125(5) 




a o 1 n nn 
4yolU.UU 


n 071 90.0.^1 fl^ 

— u.u i izyyo^iu ^ 




c;nA^9 nn 

0U4OZ.UU 


fi 7797 
— D. ii.Lt 




A7Q49 nn 
4 i y^z.uu 


ft 77/197/' 71 
— D. ( / 4Z 1 \ l ) 


9ft /I Q\ 

— zdi^ iy ) 


/iQ99ft nn 

^ioZZD.UU 


ft 771 /i 

— O. I / 14(^0 ) 




4o04o.UU 


ft 7791 z>(7\ 




aqqv^ nn 

^OO I o.uu 


-6.7716(11) 




49193.00 


-6.767 


- 


49490.00 


ft 7ft77/'Q s \ 
— O. ID / ( 1^0 ^ 




AQfizin nn 

4y04U.UU 


ft 7fi^Q/'Q s \ 

— o. * uoy \o ) 




AQS9fi nn 
^yozu.uu 


ft 7^;q/'q 1 \ 




AQQ^fi nn 


ft 771 ^n/Ql 
— Kj.it lOU^O ) 


97(V1 9^ 
Z ( U^IZ ) 


c;ni 7A nn 

QUI I 4.UU 


ft 77nnAf fi"\ 
— D. it UU4^D ) 


9A^^' r i^^ 

Z^iU^OU ^ 


c;nA9i nn 

OU4Z1 .uu 


ft 7fiQA7H s ! 


OOU^lOU ^ 


^nfin^ nn 
ouuuo.uu 


fi 771 1 ^7^1 


1 (Vl Q\ 


^n&^Q nn 
ouooy .uu 


-9.1886(3) 


-523(143) 


48409.00 


Q 1 (\^0( A\ 

— y. ioozi^ ^ 




ziQftQQ nn 

4oDOO.UU 


— y. iyuo4^4 ) 


7*3 / 1 9\ 


/1QQ9Q nn 
4oyzo.uu 


-9.1456(7) 


-211(227) 


49320.00 


-9.1552(4) 


280(180) 


49673.00 


— y.Z0044(4J 


oOU^oU J 


/Iooqi nn 
4yyoi.Uu 


O 99^09'3/'l Q\ 

— y.zzoyzoi^io ) 


— ZU4|^o j 


f^n/ioft nn 
oU4yo.uu 


-9.2199(3) 


-1600(300) 


50853.00 


Q QftQ7/IQ/'e , \ 

— o.ooo i 4y i^o ^ 


1 9zl 1 M A^i 


/ioqoo nn 

400Zo.UU 


— o.oyuui^io ^ 


101^14 ) 


4OOO0.UU 


8 &K.7AAA( ~\ 0\ 


17Q 1 ^7^ 


c;nnzi9 nn 

OUU4Z.UU 


-0.237549(7) 


- 


49204.00 


-0.237668(3) 


2.5(3) 


50059.00 


— U.ZO ( O 10(^0 ) 




r^nQ^Q nn 
ouyoy .uu 


-0.6527(6) 


- 


48357.00 


-0.653540(5) 


2.7(6) 


49054.00 


— U.0004 1 yd) 




/IQ070 nn 
4yo 1 o.uu 






^n94n nn 

OUZ4:U.UU 






^osoq on 


-8.17639(11) 


130(40) 


49141.00 


-8.19045(7) 


399(7) 


50064.00 


-8.20337(10) 


-140(60) 


50770.00 


-7.4889(3) 




50710.00 


-7.5440(6) 


4500(700) 


50831.00 


-7.52780(5) 


484(190) 


51019.00 



Span No. of Residual 



(MJD) 


TOAs 


(us) 


^nn9A ^nifiA 

OUUZ^i— OU004 


OO 


9An 

Z4iU 


50387-51155 


223 


3190 


47QnQ 4SQ9Q 

4 1 yuy— 4oyzy 


oy 


1 ni 9n 

1U1ZU 


48957-49026 


9 


113 


49043-50786 


164 


21400 


^n7Q1 ^f)QAC\ 


^1 

01 


1 on 
iyu 


ZlQI 7fi ^nzLn9 
■lyi iu ou^uz 


1 97 
iz / 


U^iZU 


0U4OO— ouoyy 


97 
Z / 


1 79 
1 1 Z 


KnAI^ ^1 1 
0U4OO— O i 1 OO 


y4 


A^nfi 

40UD 


AQ^I fi AQ7C\A 
4yolD _ 4y 1 U4 


1 ^ 
10 


QQQ 

yyo 


/iQ7ne p;i 1 
4y 1 uo _ o± 100 


1 

iyo 


si 7 

Ol / 


47Q1 ^ A7QfiQ 

4 ( y 10— 4 / yoy 


Q 


1 An9 

l^iUZ 




■^n 
ou 


1 778 
1 / / O 


ASAfi 1 ^ ASfi^l 
4o4DO — 4oDo 1 


1 fi 
10 


fi7fi 


48635-49114 


26 


2309 


49159-49227 


8 


418 


ZLQA&S AQAQ1 
4^00 — ^y^iy i 


2 




AQ^An AQ7^Q 

4yo4u— 4y 1 oy 


Q 

y 


fi78 
D 1 


AQ7fi9 AQSRQ 
^y 1 uz 'iyooa 


1 ^ 
10 


411 


AQQ9n AQQQ^ 


12 


1000 


^nn9^ ^n^9^ 

OUUZO 


^1 
01 


oyo 


^n^Ai ^n^m 

0UO41— OUOUl 


1 Q 

iy 


1 Q1 

lyi 


c;nccoft c;nfi71 

OUOOU OUU 1 -L 


9Q 
zy 


482 


50696-51022 


83 


1102 


4ozyo _ 4oOZO 


^9 
oz 


9n^A 

ZU04 


AS^Qfi ASfifiQ 

4ooyo— 4oDoy 


1 1 
1 1 


1 79 
1 ( Z 


48732-49114 


23 


820 


49159-49482 


18 


5104 


AQ^^Q AQ7S7 

4yooy— 4y / o / 


99 
ZZ 


9^AQ 
Z04y 


AQS1S A*3 
4yolo— 0U14O 


^9 
OZ 


1 nm 

1UU1 


50214-50777 


87 


1239 


50780-50926 


34 


1132 


A7QnQ A&7A& 

4 / yuy— 4o / 4u 


AA 


ooy 


ASS19 ASQ^Q 

4ooiz— ^oyoy 


91 
Zl 


loO 


Aonnn ^1 i ^fi 
4yuuu— 01 100 


9^8 
Zoo 


yoo / 


49043-49364 


24 


375 


aqai ^ ^n7n^ 
^y^io— ou 1 uo 


y^± 


^An^ 

04UO 


c;n799 ^1 1 ^fi 
OU 1 zz— OllOO 


A7 


000 


48296-48417 


11 


1451 




OO 




( ou-ouuz 1 


97 


1 1 07 
11U / 


OU11D— 0UOO4 




lUlo 


50461-51156 


100 


3246 


48896-49387 


21 


2256 


49481-50647 


89 


12262 


50656-50884 


64 


669 


50669-50751 


16 


176 


50779-50883 


13 


588 


50883-51156 


10 


152 
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Table 4. Glitch parameters for eleven pulsars 

Extrapolated Fitted 

PSR J Git. Glitch Epoch &v 3 /v A.2> s /i> Ai/ s /i/ A;> 9 /£< Q r d Resid. 

No. (Date) (MJD) (10~ 9 ) (1CT 3 ) (10" 9 ) (10~ 3 ) (d) (ps) 



0835 


-4510 


1 


961013 


50369.345(2)* 


2110(17) 


5.95(3) 


2132(52) 


0.38(2) 


916(48) 


242 


1048 


-5832 


1 


921118 


48944(2)* 


25(2) 


0.3(1) 


19(2) 






865 






2 


930216 


49034(9) 


2995(7) 


3.7(1) 


3000(10) 


0.025(3) 


100 


1630 






3 


971206 


50788(3) 


771(2) 


4.62(6) 


769(3) 


0.245(3) 


400 


465 


1105 


-6107 


1 


961130 


50417(16) 


281(3) 


1.3(2.9) 


279.7(2) 4.63(4) 






480 






2 


970613 


50610(3)* 


1.3(2) 


0.19(1) 


2.1(3) 






420 


1123 


-6259 


1 


941219 


49705.87(1)* 


749.12(12) 


1.0(4) 


749.31(14) 


0.0026(1) 


840(100) 


848 



1341 


-6220 


1 


900408 


47989(21) 


1507(1) 


0.15(6) 


1509(1) 






2352 






2 


910716 


48453(12) 


24.2(9) 


0.50(7) 


22.5(7) 


-0.51(5) 




- 






3 


920124 


486451' 1 0) 




7f1 1 

U. / (1) 






75 








4 


930527 


49134(22) 


10(2) 


0.6(2) 


13.2(13) 






- 






5 


940111 


49363(130) 


142(21) 


0.68(16) 


146(38) 


— — 


— 


— 






(5 


940620 


49523H 7 s )* 


33(3) 


— 55('9~) 


37(35) 












7 


950218 


49766(2)* 


11(11 

- LJ -V- L / 


—0.26(6) 


15(2) 












8 


950706 


49904(16) 


16(7) 


-1.9(4) 


31(1) 












9 


951018 


50008(16) 


1636(13) 


3.3(4) 


1648(3) 


0.004(1) 


300 








10 


960826 


50321.7(6)* 


27(1) 


0.61(6) 


29.9(8) 












11 


970322 


50528.9(8)* 


20(4) 


1.0(4) 


23.4(5) 












12 


970823 


50683(13) 


703(4) 


1.2(3) 


707.5(7) 








1614 


-5047 


1 


950326 


49803(16) 


6456(56) 


9.7(2) 


6460(80) 


0.538(11) 


2000 


2168 


1709- 


-4428 


1 


920605 


48778(34) 


2012.3(2) 


0.20(6) 


2028(20) 


0.133(7) 


1420(90) 


16670 


1731 


-4744 


1 


940204 


49387.3(2)* 


135(1) 


1.11(8) 


139.2(6) 


0.079(3) 


263(23) 


1171 






2 


970912 


50703(5)* 


2.6(6) 


0.8(1) 


3.1(5) 


0.25(13) 


250 




1801 


-2304 


1 


910717 


48454.1(3)* 


346(7) 


1.5(9) 


351(1) 






5061 






2 


941213 


49701(1)* 


64.7(5) 


0.18(7) 


60.8(4) 












3 


951205 


50050(5)* 


22.6(6) 


0.02(14) 


17.0(5) 












1 


961125 


50392(70) 


84(6) 


1.7(8) 


87(2) 


0.25(3) 


100 




1801 


-2451 


1 


940504 


49476(6) 


1998(7) 


4.85(28) 


1988(34) 


0.188(12) 


800(60) 


5023 






2 


970722 


50651(5) 


1237(4) 


3.87(9) 


1248(9) 


0.202(6) 


600 




1803- 


-2137 


1 


971104 


50765(15) 


3200(27) 


10.7(15 ) 


3185(25) 


0.161(6) 


855(35) 


334 






2 


971104 


50765(15) 






27(3) 


0.010(3) 


18(2) 





* Glitch epoch determined by phase fit. 



of unmodelled period irregularities affecting the timing po- 
sition. 

Extrapolation of fits to the data sets on either side of 
the glitches gives the estimates of glitch parameters listed in 
columns 5 and 6 of Table ^. Glitch parameters from TOA 
fits are given in the remaining columns. Because of the con- 
taminating effect of the period fluctuations, the parameters 
of the large glitches were determined by fitting the data span 
100 - 150 days before and after each glitch. The effects of 
the earlier small glitch were determined separately by fitting 
the interval between 80 days prior to it and up to the large 
glitch, and were subtracted from the fit to the large glitch. 
Parameters from these fits are given in Table Because of 
the large systematic period variations, it is not possible to 
reliably measure the post-glitch decay times. Fig. [] gives 
some evidence for a relaxation with a timescale of ~ 100 d 



following the 1993 glitch and somewhat longer for the 1997 
glitch. Setting the decay time to 100 d and 400 d respec- 
tively for these two glitches gave the parameters listed in 
Table ^[ The derived Q for the 1997 glitch appears signifi- 
cantly larger than that for the 1993 glitch, but this result is 
uncertain given the uncertainty in the relaxation time and 
the short data span following the glitch. 



3.3 PSR J1105 6107 

PSR J1105— 6107 is a young pulsar with the relatively short 
period of 63 ms discovered in 1993 using the Parkes radio 
telescope ( Kaspi et al. 1997 ). It is located close to but out- 
side the boundaries of the supernova remnant G290. 1—0.8 
and within the erro r box for the EGR ET gamma-ray source 
2EG J1103-6106 ( |Kaspi et al. 2000|) . Its association with 
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Figure 2. Glitches of PSR J1048-5832: variations of (a) fre- 
quency residual Av relative to the pre-glitch solution, (b) an ex- 
panded plot of Av where, for the second and third glitches, the 
mean residual up to the next glitch, or to the end of the data 
if there is no following glitch, has been removed from the corre- 
sponding interval (indicated by the raised arrows), and (c) the 
variations of v. 
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Figure 4. Glitches of PSR J1105-6107: variations of (a) fre- 
quency residual Av relative to the pre-glitch solution, (b) an ex- 
panded plot of Av where the mean residual between the large 
glitch and the second glitch has been removed from the data fol- 
lowing the large glitch (indicated by the raised arrow), and (c), 
the variations of v. 
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Figure 3. Timing residuals for PSR J1048-5832 between the 
two large glitches. 



these objects remains plausible but unproven. Kaspi et al. 
(1997) reported on timing observations from 1993 July to 
1996 July. Here we extend this interval to the end of 1998, 
and show that there were two glitches, a large one around 
the end of November, 1996 (MJD ~ 50417) , and a small 
one about 200 days later. 

Fig. ^ gives the variations in pulse frequency and fre- 
quency derivative over the observed data span, clearly show- 
ing the larger glitch which is of fractional size Av g /v ~ 
0.28 x 10~ 6 . There was a gap between bracketing observa- 
tions of about 30 days, so the glitch epoch is not well de- 
termined. The expanded plot (b) shows that, despite the 
substantial period irregularities, there is good evidence for 
a small glitch about 200 days after the large glitch. Fig. ^ 
shows phase residuals around the time of this small glitch. 
Although this glitch is by far the smallest discussed in this 
paper, there is little doubt about its reality. Another small 
glitch may have occurred near the end of the data set, but 
there are insufficient post-glitch observations to confirm this. 
The position given in Table ^ was determined from a fit 
to the data prior to the large glitch, with seven frequency 
derivatives to absorb the period irregularities. 

Table ^ gives fits to the three data spans delimited by 
the data set and the two glitches; the position was held at 
the value quoted in Table |^. The second derivative terms 
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Figure 5. Timing residuals around the time of the small glitch 
in the period of PSR J1105-6107. 



in the first and third fits are dominated by random period 
irregularities. Glitch parameters obtained by extrapolating 
these fits and from phase fits across the glitches are given 
in Table Because of the period noise, the phase fits were 
restricted to intervals of 150 days before the first glitch, the 
inter-glitch interval, and 150 days after the second glitch. 
The fractional increase v after the first glitch is small, im- 
plying a similarly small value of the ratio Q/tj (Equation^). 
Because of the small data span, t c i is not well determined, 
but taking 100 d gives Q ~ 0.035; i.e., only a few percent of 
the glitch is likely to decay. This is confirmed by a fit to the 
whole data span including both glitches. 
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Figure 6. The 1994 December glitch of PSR J1123-6259: vari- 
ations of (a) frequency residual Au relative to the pre-glitch so- 
lution, (b) an expanded plot of Ai> where the mean post-glitch 
value has been subtracted, and (c) the variations of i>. 



3.4 PSR J1123 6259 

PSR J1 123— 6259, discovered in the Parkes southern pulsar 
survey (Manchester et al. 1996), has a period of 271 ms and 



a modest period derivative, giving a characteristic age of 
~ 8x 10 5 yr. Despite this relatively large characteristic age, a 
glitch was detected between 1994 December 17 - 22 (MJD ~ 
49706). The pulsar position was determined from the post- 
glitch data (Table ^); no pre- whitening was necessary as 
random period irregularities are small in this pulsar. Fits to 
the pre- and post-glitch data with this position are given 
in Table ^. Extrapolation of these fits to a central epoch 
(Table [lj) shows that the glitch was of magnitude Avjv ~ 
7.5 x 10~ 7 . | 

As Table ^ and Fig. ^ show, there was a small but sig- 
nificant increase in \v\ at the time of the glitch. A fit to 
the entire data span including glitch parameters but with 
the position and pre-glitch frequency parameters held to the 
values given in Table ^ and Table ^| is an excellent fit to the 
data and gives an unambiguous value for the epoch of the 



3 A glitch of this magnitude was reported by D'Amico et al. 
(1998) to have occurred at MJD 48650 ± 20 (1992 December). 
This epoch is in error. It was in fact the December 1994 glitch 
discussed in this paper. 



glitch with an uncertainty of about 15 min on the assump- 
tion that there was phase continuity across the glitch. The 
derived value of Q is very small and the decay time is long 
(Table 0). 



3.5 PSR J1341 6220 

PSR J1341— 6220 is a young pulsar (characteristic age ~ 
12,000 y) discovered by Manchester et al. (1985) and asso- 
ciated with the supernova remnant G308— 0.1 by Kaspi et 
al. (1992). Kaspi et al. reported two large glitches (Av g /u ~ 
1.5 x 10" 6 and ~ 1.0 x 10~ 6 ) and one smaller one (~ 
2.3 x 10~ 8 ) within the period 1990 January to 1992 May, 
making this one of the most actively glitching pulsars known. 
In this paper, we reanalyse the data from 1990 to 1992 and 
extend the data set 1998 March. Unfortunately, observations 
after this time were too sparse to permit unambiguous pulse 
counting. 

The position given in Table ^ was determined in a simul- 
taneous fit across the whole data set, including all glitches 
(see below). It has smaller estimated uncertainties than the 
position given by Kaspi et al. (1992), and lies 4'.'3 north of it. 
Given the prominence of period irregularities in this object, 
this difference, while twice the combined uncertainties, is of 
marginal significance. Determination of the DM is compli- 
cated by the strong scattering tail shown by this pulsar, even 
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Figure 7. The glitches of PSR J1341-6220: variations of (a) 
frequency residual Au relative to the pre-glitch solution, (b) an 
expanded plot of Av where the mean residual between glitches 
with a raised arrow and the following glitch has been removed 
from data after the marked glitch, and (c) the variations of v. 

at frequencies around 1400 MHz. The value given in Table ^ 
was determined using TEMPO with data recorded between 
1998 March and December at frequencies close to 1400 MHz 
and 1700 MHz. The relative phase of the standard profiles 
at the two frequencies was adjusted to allow for the effects 
of scattering. The derived DM is much more accurate than 
the value given by Kaspi et al. (1992) and just outside the 
error range of this value. 

Fig. M shows the very complicated period history of this 
pulsar. A total of 12 glitches were observed in the 8.2-year 
interval, making this the most frequently glitching pulsar 
known, with a mean glitch interval of 250 days. The next 
most freq uent is PSR B1737— 3 0, with a mean interval of 
345 days (Shemar & Lyne 199£). Four of the glitches were 
relatively large, with Av g /v > 0.7 x 10~ 6 and for two the 
fractional glitch size exceeds 10~ 6 . Most of the other glitches 
are quite small, with all except one having Au g /u < 0.04 x 
10 -6 . They are, however, unambiguous on phase-time plots. 

For most of the glitches, there is little evidence for any 
post-glitch relaxation. This is partly because of the relatively 
poorly sampled data, especially in the earlier years, but also 
because of the short time between successive glitches. For 
the third and ninth glitches, there is weak evidence in the i> 
plot for some relaxation. 

Table ^ gives fits to the inter-glitch intervals. Because of 
the limited data span prior to the first observed glitch, the 



value of i> was held fixed at a value close to the long-term 
mean for this fit. For about half of the intervals, there was 
a clear v term in the residuals and this term was fitted. 

Changes in v and v at the time of each glitch, obtained 
by extrapolation of these fits to the glitch epoch, are given 
in Table ^. This table also gives glitch parameters from a 
single fit to the whole data set, solving simultaneously for 
the pulsar position, pulsar frequency (u), mean frequency 
derivative (i>), and the parameters for all 12 glitches. As 
mentioned above, for most of the glitches, there was no 
significant post-glitch relaxation. For the third and ninth 
glitches, a post-glitch exponential decay was fitted; in both 
cases the decay time constant was not fitted for, but was 
determined by trial to minimise the final residual. In both 
cases the derived Q values are relatively small. 

3.6 PSR J1614 5047, PSR B1610 50 

This pulsar has a period of 231 ms and a very large period 
derivative (495 x 10 -15 ), implying a small characteristic age 
of ~ 7400 yr. In terms of its period irregularities, it is one of 
the noisiest pulsars known; this makes it difficult at times to 
keep track of the puls e phase. Since its discovery in late 1989 
(Johnston et al. 1992), there has been no clear evidence for a 
glitch although, particularly where there is a significant gap 
in the timing data, it is often difficult to distinguish between 
a (small) glitch and more continuous period irregularities. 
However, in 1995 June (MJD ~ 49802), there was a massive 
glitch with Avg/v ~ 6.5 x 10~ 6 , the largest ever observed in 
any pulsar (cf. Shemar & Lyne 1996). Fig. |^ shows this glitch 
and also illustrates the more continuous period irregularities. 

These irregularities make it difficult to determine the 
pulsar position from timing data. Previously published po- 
sitions h ave differed by much more than the quoted uncer- 



taint ies (Taylor, Manchester fc Lyne 1993 Johnston et al 



199J). Taking the data span from MJD 50269 to to 50778 
(which is free of major irregularities - see below) and fitting 
for position and two frequency derivatives gives the position 
quoted in Table ^. An independent fit to the MJD range 
48732 - 49093 gave a position with uncertainty of about 2", 
consistent with the position from the longer data span. Sub- 
sequent fits were made keeping the position fixed at the Ta- 
ble ^ value. Stappers et al. (1999) have recently determined 
an interferometric position for this pulsar: R.A. (J2000) 16 h 
14 m 11:55=^01, Dec. (J2000) -50° 48' 01'.'9 ± Cf.'l. As for 
PSR J1048— 5832, the differences most probably result from 
period irregularities affecting the timing position. 

The slope of the post-glitch data in the top two plots in- 
dicates a large change in frequency derivative at the time of 
the glitch. Because of the strong irregularities in the period, 
the pre- and post-glitch fits given in Table |^ are restricted 
to intervals of less than 300 days. Extrapolation of those fits 
gives the glitch parameters in columns 5 and 6 of Table ^j. 
Glitch parameters were also obtained from a single timing 
solution over the same interval as the pre- and post-glitch 
fits (MJD 49559 - 50117) and are given in columns 7 - 10 of 
Table ^. This fit did not converge well when fitting for the 
glitch decay time, so this was held fixed at 2000 d, a value 
representative of those obtained from the fitting. From the 
post-fit residuals, it is clear that there is also a more rapid 
decay of a portion of the glitch with a timescale of 10 - 20 
days. There is good agreement between the two methods of 



10 



PSR J1614-5047 



o 

K 

o 

<l 



a. 



< 



M O 



0) 



1 1 

: (a) 






• : 




• • • • • • 


• « 




— 1 — 

- (b) 


H 1 1 1 


- 

• 

• 


H 1 1 1 ; 


• » % 




• • 

T 


* 

% — 
• 


— 1 — 

" (c) 

— V / 
m 

m 


H 1 1 ( 

+ 

• • 


1 1 — 

• 


H 1 1 1 ; 


— • 

1 




• 

t 


' . .. 





1000 




2000 



MJD 



48300 



Figure 8. The 1995 June glitch of PSR J1614-5047: variations 
of (a) frequency residual Au relative to the pre-glitch solution, 
(b) an expanded plot of Av where the mean residual after the 
glitch has been removed, and (c) the variations of v. 



deriving the glitch parameters. From Equation ^, the Q and 
Td values given in Table ^ imply a fractional change in v of 
0.0096, compared to 0.0097 from the extrapolation. Unfor- 
tunately, because of the large gap between the observations 
bracketing the glitch and the large size of the glitch, it is 
not possible to accurately determine the glitch epoch from 
phase continuity. 

This is the largest glitch yet observed. The previ- 
ous largest w as for PSR B0355+54 for which Av/v was 
~ 4.37x10" 



(Lyne 1987). There was very little decay of the 
frequency step in PSR B0355+54, whereas the Q value for 
the PSR J1614-5047 glitch indicates that more than half 
of it will decay away over a few years. In this respect, the 
PSR J1614— 5047 glitch is more similar to the giant glitches 
in the Vela pulsar, which decay by a similar amount over 
similiar timescales. 

Phase fitting of post-glitch data revealed two further 
abrupt changes in pulse frequency. Phase plots for the inter- 
vals around these events are shown in Fig. ^. These phase 
plots have the character of glitch events, with a persistent 
fractional frequency change Av g /v ~ 0.03 x 10 -6 at about 
MJD 50170 and 50780, respectively. The first event may in- 
deed be a small glitch - it is not possible to tell because of 
the large data gap - but the second is not. This frequency 
change is resolved, taking place over about 10 days. In both 
cases, there appears to be a significant decrease in \v\ asso- 



ciated with the event. These rapid frequency changes could 
be classed as just part of the frequency irregularities which 
are prominent in this pulsar, but in general these irregular- 
ities have a much longer timescale and are smaller in am- 
plitude. This is demonstrated by the small residuals for the 
fits to data before these events; for the second event, the 
fit includes only two frequency derivatives and has an rms 
residual of only 1.2 ms (Table |). The MJD 50780 event 
stands out from the general irregularities with a much larger 
rate of frequency change over the 10 days. Phase fits to the 
post-glitch data given in Table [] have been split into three 
sections to avoid contamination by these sudden frequency 
changes. The last ends at MJD 50926 (1998 April) since data 
are sparse after that point and it is not possible to fit across 
the gaps without ambiguity. The frequency second deriva- 
tive terms are very significant in these fits, but differ greatly 
in both sign and magnitude for the different segments. They 
are clearly related to the period-noise processes occurring in 
this star and not to the secular slowdown. 



3.7 PSR J1709 4428, PSR B1706-44 

This pulsar, discovered by Johnston et al. (1992), is of in- 
terest for several reasons. It is young (r c ~ 17 kyr) and 
has the third highest known value of E/d 2 , where E is the 
spin-down luminosity and d is the pulsar distance, after the 
Crab and Vel a pulsars. It has been detected at gamma-ra y 
wavelengths ( Thompson et al. 1992 ; Thompson et al. 1996 ), 
X-ray wavelengths (|Be'cker7^razier^r*Triu*nper 1995[ ) and 
possibly at TeV energies ( |Kifune et al. 19950 . A' supernova 
remnant association was proposed by McAdam, Osborne & 
Parkinson (1993), but deemed unlikely by Nicastro, John- 
ston & Koribalski (1996). 

Johnston et al. (1995) presented results from two years 
of Parkes timing observations showing that the pulsar suf- 
fered a giant glitch of magnitude Av g /u ~2x 10~ 6 near the 
end of May, 1992 (~ MJD 48778). In this paper we analyse 
data from 1990 January to 1998 December. To determine 
the pulsar position, we have taken data from 150 days after 
the glitch (MJD 48928 - 51155) and fitted nine frequency 
derivatives to absorb the period irregularities. This fit had 
an rms residual of only 211 /xs and gave the position listed 
in Table |i| The derived position is about 12" southeast of 
the position given by Johnston et al. (1995), but agrees bet- 
ter with an interferometric position given by Frail & Schar- 
ringhausen (1997): R.A. (J2000) 17 h 09 m 4275 ± K01, Dec. 
(J2000) -44° 29' 06'.'6±0'.'3. The timing fits described below 
were obtained with the position held at the Table ^ value. 

Fig. [H^ shows the large jump with a classic exponential 
recovery on a timescale of ~ 100 days coupled with a longer 
term relaxation. Fits of a cubic phase polynomial to the pre- 
glitch data, to the 150 days after the glitch and to data from 
150 days post-glitch are given in Table [| The post-glitch fit 
has large systematic residuals, dominated by the post-glitch 
recovery. Glitch parameters obtained by extrapolating the 
first and second fits to the glitch epoch are given in Table |j| 
along with those from a phase fit to the whole data set. 
This latter fit is strongly affected by the strong period ir- 
regularities (Fig. [ll]) and so the derived parameters are only 
approximate. There is clear evidence in the post-fit residu- 
als for a more rapid exponential decay with time constant 
of the order of 100 days. However, it was not possible to fit 
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Figure 9. Variations of timing phase residuals around two epochs of rapid frequency change in PSR J1614— 5047 
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Figure 10. The 1992 May glitch (MJD 48778) of PSR 
J1709— 4428: variations of (a) frequency residual Av relative to 
the pre-glitch solution, (b) an expanded plot of Av where the 
mean residual after the glitch has been removed, and (c) the vari- 
ations of v. 



for the parameters of this owing to the contaminating effect 
of the period irregularities. 

3.8 PSR J1731-4744, PSR B1727-47 

PSR J1731-4744 is a long-period pulsar (0.830 s) which has 
a large period derivative giving a relatively low characteristic 
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Figure 11. Post-fit timing residuals for a fit to the entire data 
set for PSR J1709-4428, including the glitch at MJD 48778 



age, ~ 80,000 years. D'Alessandro & McCulloch (1997) ob- 
served this pulsar at Mt Pleasant Observatory in Tasmania 
from 1987 to 1994 and detected a relatively large glitch at 
MJD 49387.68 (1994 February 2) but found no evidence for 
any recovery in 10 months of observations after the glitch. 
The Parkes observations extend to the end of 1998 (Fig. |l2[ ) 
and show a clear recovery from this glitch which is roughly 
exponential. They also reveal a second smaller glitch in 1997 
September which has a similar recovery. 

Table ^ gives an improved position and DM for this 
pulsar. The position was determined from data between the 
two glitches by fitting for the position and five frequency- 
derivative terms to absorb the period irregularities. The DM 
was determined by fitting to data in the MJD range 49400 
to 49900 where there were observations around 430 MHz as 
well as around 1400 MHz. This value is significantly differ- 
ent from the best previously published value (121.9 ± 0.1 
cm -3 pc; McCulloch et al. 1973), most probably reflecting a 
changing line-of-sight path through the interstellar medium. 
The average rate of DM change over the 24 years is ~ 0.6 
cm -3 pc y~ . DM changes have been observed in other pul- 
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Figure 12. Glitches of PSR J1731-4744: variations of (a) fre- 
quency residual Av relative to the pre-glitch solution, (b) an ex- 
panded plot of Av where the mean residual in the inter-glitch 
interval has been removed from the following data, and (c) the 
variations of v. 



s ars, notably the Crab ( [ tsaacmaii fc Rankin 1977 ) and Vela 



(Hamilton, Hall & Costa 1985) pulsars. The average rate of 
change of DM for PSR J1731-4744 is comparable to that 
observed for the Vela pulsar and larger than that for the 
Crab pulsar. The large changes observed for the these two 
pulsars are attributed to ionised gas within the associated 
supernova remnant, but PSR J1731— 4744 has no associated 
supernova remnant. 

Fits to the inter-glitch intervals are given in Table |!| 
These fits represent the data reasonably well, except that 
there is a significant quartic term in the residuals for the 
middle interval. This is mostly due to the period irregulari- 
ties which are evident in Fig. [l2| 

Table ^ gives parameters for the two glitches. The fitted 
parameters were determined by a simultaneous fit to both 
glitches. As mentioned above, a recovery after both glitches 
is clearly seen in Fig. [L2| and these are well fitted by the ex- 
ponential model. For the first glitch, about 8 per cent of the 
glitch is recovered with a time constant of about 260 days. 
This implies a fractional increment in v at the time of the 
glitch of 2.25 x 10 -3 . The glitch parameters are in reasonable 
agreement with those quoted by D'Alessandro & McCulloch 
(1997); the differences of a few times the combined uncer- 
tainties could result from the different models used for the 
extrapolation to the time of the glitch. The data span fol- 
lowing the second glitch is too short to permit solving for the 



time constant, so 250 days was assumed, giving a Q value 
of about 25 per cent. 



3.9 PSR J1801 2304, PSR B1758-23 

This pulsar was detected in a search f or short-period pulsars 
associated with supernova remnants (Manchester, D'Amico 
& Tuohy 1985). Despite its relatively long period of 415 

It 



ms, this pulsar has a short characteristic age (58 kyr). 
lies close to the supernova remnant W28, but its association 



with the remnant remains cont roversial (Kaspi et al. 199S 



Frail, Kulkarni & Vasisht 1993). It has the highest known 
dispersion measure (1074 cm~ d pc) and even at 1.4 GHz the 
profile is very scattered, reducing the precision of pulse tim- 
ing observations at this and lower frequencies. Furthermore, 
the pulsar lies very close to the ecliptic plane and suffers 
frequent glitches, so that the pulsar position is not well de- 
termined by pulse timing; the position quoted in Table |^ and 
used for the timing analyses is from the VLA observations 
of Frail et al. (1993). 

Four glitches in this pulsar occurring between 1986 and 
the end of 1994 were reported by Kaspi et al. (1993) and 
Shemar & Lyne (1996). In this paper, we present data on 
the last two of these glitches and two further glitches, both 
of relatively small amplitude, occurring in 1995 and 1996 
respectively. Fig. shows the variations in pulse frequency 
and frequency derivative around these four glitches relative 
to the pre-glitch solution given Shemar & Lyne (1996) for 
the MJD 48454 glitch. Thispre-glitch solution was adopted 
because that given in Table tfl is based on a data span of only 
120 days and has a significantly different value of i> from the 
other fits, probably as a result of intrinsic period irregular- 
ities. Fits to the inter-glitch intervals given in Table N have 
relatively large residuals because of these irregularities. Only 
the fit to the data following the fourth glitch required a v 
term, indicating a significant relaxation following the glitch. 

Parameters for the four glitches are given in Table ^j, 
both from extrapolation of the polynomial fits given in Ta- 
ble ^ and from a simultaneous fit to the TOA data across 
all four glitches. Values of Av/v and Avji> for the first two 
jumps have smaller uncertainties but are consistent with the 
values given by Shemar & Lyne (1992). For the first three 
jumps the glitch epoch is determined from the requirement 
of phase continuity over the glitch; for the fourth glitch the 
data gap is too large for an unambiguous determination of 
the glitch epoch. Consistent with the observation of a sig- 
nificant v term, post-fit residuals were significantly reduced 
by allowing an exponential relaxation following the fourth 
jump. However, the presence of period irregularities made 
fitting for the decay time impossible; the value of 100 days 
was found by trial to give the best representation of the 
data. 



3.10 PSR J1801 2451, PSR B1757-24 

PSR J1801-2451 is a young pulsar (r c ~ 15, 000 y) which is 
located at the "beak of the Duck", in the small nebula as- 
s ociated with the much l a rger supernova remnant G5.4— 1.2 



(Manchester et al. 1991; Frail fc Kulkarni 1991). A giant 



glitch in the pulse period was observed around MJD 49476 
by Lyne et al. (1996). About a year of post-glitch data 
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Figure 13. Glitches of PSR J1801-2304: variations of (a) fre- 
quency residual Av relative to the pre-glitch solution, (b) an ex- 
panded plot of Av where the mean residual following each glitch 
with a raised arrow has been removed, and (c) the variations of 
v. 
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Figure 14. Glitches of PSR J1801-2451: variations of (a) fre- 
quency residual Av relative to the pre-glitch solution, (b) an ex- 
panded plot of Av where the mean residual after each glitch has 
been removed, and (c) the variations of v. 



to MJD 49950 (1995 August) suggested a relaxation of a 
small fraction of the glitch (Q ~ 0.005) with a characteristic 
timescale of about 40 days. 

In Fig. [ll| we show observed pulse frequency variations 
based on Parkes data from 1992 October to 1998 Decem- 
ber. These show the glitch observed by Lyne et al. (1996) 
and also a second smaller, but still large, glitch around 1997 
July. Polynomial fits to the inter-glitch intervals given in 
Table ^| have large residuals owing to the presence of period 
irregularities. Extrapolation of the first two of these fits to 
the glitch epoch determined by Lyne et al. (1996) gives re- 
sults consistent with those obtained by these authors, with 
Avg/u ~ 2.0 x 10~ 6 . The second glitch, separated from 
the first by about 3.2 years, has a fractional amplitude of 
~ 1.2 x 10 -6 , making it also a member of the giant glitch 
class. 

The longer timespan of the present observations follow- 
ing the first glitch show that the post-glitch relaxation is 
best described by an exponential recovery with a timescale 
of several hundred days. It is likely that period irregularities 
were primarily responsible for the apparent quicker relax- 
ation found by Lyne et al. (1996). There is also evidence 
for a similar long-timescale relaxation following the second 
glitch. A TEMPO fit across both glitches yielded the pa- 
rameters given in Table ^, showing that, for both glitches, 
about 20 per cent of the glitch decayed. For the first glitch, 



the fitted timescale was 800 days. No fit of the timescale to 
the second glitch was possible because of the shorter data 
span; an assumed decay time of 600 days gave a minimum 
in the post-fit residuals. 



3.11 PSR J1803 2137, PSR B1800 21 

PSR J1803-2137 is a young pulsar (r c ~ 16,000 y) dis- 
covered by Clifton & Lyne (1986). Lyne & Shemar (1996) 
observed a large glitch (Av g /v ~ 4.1 x 10 -6 ) in 1990 Decem- 
ber which showed an exponential recovery with timescale of 
~ 150 d, together with an apparently linear decay in v in- 
dicating decay from an earlier unobserved glitch. 

Our observations of this pulsar are from 1997 August 
to 1998 December. Fig. ^ shows that another large glitch of 
magnitude Ai/ g /v ~ 3.2 x 10 -6 occurred in 1997 November. 
Again, there the signature of an exponential decay following 
the glitch is seen. There is, however, no evidence of decay 
in z> preceding the glitch (Table although the data span 
is rather short. Table ^ gives the extrapolated and fitted 
parameters for the glitch. Fitting of a single exponential de- 
cay gives a relatively good fit to the post-glitch data with 
Q ~ 0.13 and Td ~ 640 d and an rms residual of 1690 /is. 
However, there is clear evidence in the residuals immedi- 
ately after the glitch for a shorter-term decay. The glitch 
parameters in Table U are from a simultaneous fit of two ex- 
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Figure 15. The 1997 December glitch of PSR J1803-2137: vari- 
ations of (a) frequency residual Av relative to the pre-glitch so- 
lution, (b) an expanded plot of Av where the mean residual after 
the glitch has been removed, and (c) the variations of v. 



ponential decays, one with a short decay time and the other 
representing the longer-term decay. The final fit is extremely 
good, with a final rms residual of only 334 fis and the resid- 
uals dominated by random noise. Fitting of the short-term 
decay resulted in an increase in the estimates for Q and Td 
for the longer-term decay by about 25 per cent as shown in 
Table §. 



4 DISCUSSION 

In this paper we have presented timing observations of 
40 pulsars over various intervals up to a maximum of 8.9 
years and analysed these data for improved astrometric and 
pulse frequency parameters and for glitch activity. In to- 
tal, 30 glitches were detected in 11 pulsars, including the 
largest known glitch, with Av g /v ~ 6.5 x 10~ 6 in PSR 
J1614— 5047. Twelve glitches were detected in the period 
of PSR J1341— 6220 in a data span of 8.2 years, making this 
the most frequently glitching pulsar known. Evidence was 
found in PSR J1614— 5047 for a new class of irregularity in 
which the pulse frequency increases markedly over a few-day 
interval. There appears to be an accompanying decrease in 
the magnitude of the frequency derivative, implying a cor- 
respondingdecrease in braking torque. 

Table n lists all known glitches, 76 in total, giving the 
fractional glitch amplitude Av g jv and, in parentheses, the 



approximate MJD of the glitch. Fig. |l6| shows a histogram 
of fractional glitch amplitudes. Largely because of the Vela 
pulsar, the most common glitches are large, with Av g /v in 
the range 1 — 3x 10~ 6 . In most other pulsars, smaller glitches 
with fractional amplitude down to 10 -8 are more common. 
Glitches with size smaller than 10 -9 are difficult to identify, 
especially in noisy pulsars, and the sample is certainly in- 
complete at this level. However, there does appear to be a 
reduced rate of occurence of glitches with Au g /u < 10 -8 . 

In most models for the glitch phenomenon, the sudden 
spin-up is triggered by the release of stress built up as a 
result of the stead y spin down of the pulsar. For the original 
star-quake model (Baym et al. 1969), the equilibrium shape 
of the star becomes less oblate as the star spins down. At 
some point, the crust cracks and relaxes to (or toward) the 
new equilibrium shape, reducing its moment of inertia and 
hence spinning it up. However, this model fails to account 
for frequent giant glitches, as seen, for example, in the Vela 
pulsar, as the rate of change of oblateness is too slow. In 
models based on unpinning of internal superfluid vortice s 
flAlpar, Cheng fc Pines 1989; |Ruderman, Zhu fc Chen 1998[), 



the stresses on the pinned vortices build up as the crust slows 
down until finally some fraction of them unpin and then 
repin at a larger radius, resulting in a transfer of angular 
momentum to the crust. 

In these models, one expects some relation between the 
size of the glitch and the length of time that the pulsar 
has been slowing down since the previous glitch or the in- 
tegrated change in spin rate since the last glitch. Fig. |l7|(a) 
shows glitch fractional amplitudes plotted against the length 
of the preceding inter-glitch interval (Ai 9 ), and in Fig. p"7|(b) 
against the change in spin frequency since the previous glitch 
(|i»|Aig). These figures show that, contrary to expectations, 
there is no general relation between fractional glitch ampli- 
tude and either the time since the previous glitch or the 
total change in spin frequency since the previous glitch. For 
the Vela pulsar (marked with a * in the Figure), the smaller 
glitches do tend to have shorter preceding intervals and all 
but one of the giant glitches have preceding inter-glitch inter- 
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Table 5. Summary of all known glitches 
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vals of about 1000 days. However, for PSR J1341-6220 there 
is if anything an inverse relationship, with larger glitches oc- 
curing after shorter intervals. For PSR J1740— 3015 there is 
no relationship between glitch size and preceding interval. 
Fig. 0(b) shows that similar relationships for individual 
pulsars hold when glitch size is plotted against accumulated 
spin-down since the last glitch. The three points to the right 
on this plot are for the Crab pulsar; these too show an in- 
verse relationship. 

Alternatively, if a small glitch results from a release of 
only part of the built-up stress, one might expect another 
glitch to occur soon after, when the breaking strain is again 
reached. Conversely, if a large glitch releases all or most of 
the stress, a long time would be required for it to build up 
again. This suggests a correlation between glitch size and 
time to the next glitch. Fig. |l^ shows glitch size plotted 
against duration of the following inter-glitch interval and 
accumulated spin-down in that period. No such correlation 
is observed for the Vela pulsar, where the small glitches are 
followed by relatively long intervals, but a weak positive cor- 
relation is seen for PSR J1341-6220 and PSR J1740-3015. 
Excepting the Crab pulsar, other pulsars (marked with a 
dot) have a good correlation between glitch size and accu- 
mulated spin-down frequency following the glitch. 

These results suggest that the triggering of glitches is 



a local phenomenon, not dependent on global stresses. This 
lends some support to the ideas of Ruderman et al. (1998) 
in which migration of magnetic flux tubes determines the 
stresses on the pinned vortices. This model is also supported 
by the observations of the "slow" glitches and subsequent 
decrease in slow-down rate in PSR J1614— 5047. 

McKenna & Lyne (1990) introduced the glitch activ- 
ity parameter, A g , defined to be the accumulated pulse fre- 
quency change Av g due to glitches divided by the obser- 
vation data span. It therefore has the same units as v and 
represents the portion of v which is overcome by glitches. 
This is typically small; the largest known value of A g , for 
the Vela pulsar, is ~ 2.5 x 10 -13 s -2 , only about 2 per cent 
of its spin-down rate v ~ 1.6 x 10 -11 s~ 2 . Fig. [l9|(a) shows 
activity parameter plotted against characteristic age, t c for 
most of the pulsars with extensive timing data. As noted by 
McKenna & Lyne (1990), there is a clear peak in activity 
for pulsars with ages between 2,000 and 20,000 years. 

However, there is a substantial group of young pulsars 
with low glitch activity; these so far have had no observed 
glitch and are represented by upper limits in Fig. [19(a) cor- 
responding to a single glitch of fractional size 10 . These 
eight pulsars have a mean data span of more than six years 
and so, at the least, they are infrequent glitchers. A single 
glitch of amplitude ~ 10 -6 would raise them into the same 
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Figure 17. (a) Fractional glitch amplitude Au g /u versus length of the preceding inter-glitch interval, (b) Fractional glitch amplitude 
Av g /v versus accumulated spin-down frequency during the preceding inter-glitch interval. In both plots, points for the Vela pulsar (PSR 
J0835-4510) are marked with a *, those for PSR J1341-6220 are marked with +, and those for PSR J1740-3015 (PSR B1737-30) are 
marked with o. 
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Figure 18. (a) Fractional glitch amplitude Au g /u versus length of the following inter-glitch interval, (b) Fractional glitch amplitude 
Av g /v versus accumulated spin-down frequency during the following inter-glitch interval. In both plots, points for the Vela pulsar (PSR 
J0835-4510) are marked with a *, those for PSR J1341-6220 are marked with +, and those for PSR J1740-3015 (PSR B1737-30) are 
marked with o. 



region as the other young pulsars. The dashed line corre- 
sponds to a uniform accumulated v over the characteristic 
spin-down time. This is approximately true for pulsars with 
age of less than or about 10 6 years, but older pulsars clearly 
have less glitch activity than predicted by this rule. 

In Fig. |l^(b) the activity parameter is plotted against 
the absolute value of v. Some theoretical models for glitches 
(e.g. Ruderman, Zhu & Chen, 1998) predict that the ratio 
of the effective spin-up rate due to glitches to the spin-down 
rate is proportional to the ratio of the moment of inertia 
of the crustal superfluid to the total moment of inertia of 
the neutron star. A constant ratio corresponds to a line 
of slope —1 in Fig. ^(b). For the younger pulsars, except 
PSR B0531+21 (Crab) and PSR B1509-58, the points cor- 



respond to a fraction of 1 - 2 per cent (cf. Lyne et al. 1999). 
Pulsars with spin-down rates less than about 10~ 14 s _1 have 
a significantly smaller fraction of their spindown recovered 
by glitches. 

When a significant post-glitch decay is observed, it is 
generally well described by the exponential relation given in 
Equation |[ In Fig. |2o|(a) the fraction of the glitch which 
decays, Q — Aval ' tAv g is plotted against characteristic age. 
With one exception, PSR B0525+21, old pulsars tend to 
have low values of Q. Younger pulsars can have larger Q, 
but many glitches in young pulsars have little or no decay. 

Fig. M(b) is a plot of glitch decay timescale versus pul- 
sar characteristic age. This shows that, apart from the Crab 
pulsar (the three points in the lower left of the figure) in 
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whereas those with greater r c > 10 5 yr are mean values in half-decade bins. Where no glitch was observed, the upper limit is based on a 
single glitch of fractional amplitude 10~ 9 over the total data span, summed over all pulsars in the case of binned data. Upper limits are 
higher for the very old pulsars simply because there are fewer of them known and the total accumulated data span is less. The dashed 
line has a slope of —1.0. (b) A similar plot with the absolute value of the spindown rate i> as abscissa. Points with \v\ > 10~ 12 are plotted 
separately while those with \v\ less than this value are binned. 
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Figure 20. (a) Plot of Q, the fraction of Ai/ g which decays, versus characteristic age. (b) Plot of glitch decay timescale versus characteristic 
age. Results from this work are marked with a *. 



which the decay timescale is very short, there is no relation 
between decay timescale and pulsar age. We emphasise that 
this result applies to the longest decay timescale present. 
This may indicate the unreliability of characteristic age as 
an indicator of true age and hence internal temperature, or it 
may indicate that decay time, at least for pulsars with age 
greater than a few thousand years, is dependent on other 
properties, for example, the core magnetic field. 

These observations have demonstrated the great diver- 
sity of glitch properties. Glitch activity is clearly greatest 
in pulsars with ages of between a few times 10 and 10 
years. The Crab pulsar has distinctly different glitch prop- 
erties from those of the middle-aged pulsars, and other pul- 



sars of similar age show no glitches at all. Apart from these 
clear trends, glitch properties vary greatly, both for succes- 
sive glitches from a given pulsar and across different pulsars, 
with few systematic trends. These properties suggest that 
the glitch phenomenon, both the event and the following re- 
sponse, depend on quasi-random processes occurring in the 
pulsar, rather than global properties such as slow-down rate 
or age. Crustal processes driven by magnetic field evolution 
as proposed by Ruderman et al. (1998) seem more consis- 
tent with this than the vortex creep models of Alpar et al. 
(1989). 
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